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ABSTRACT 

We present reduced and calibrated high resolution Lyman-alpha (Lya) images for a sample 
of six local star forming galaxies. Targets were selected to represent a range in luminosity and 
mctallicity and to include both known Lya emitters and non-emitters. Far ultraviolet imaging 
was carried out with the Solar Blind Channel of the Advanced Camera for Surveys on Hubble 
Space Telescope in the F122M (Lya on-line) and F140LP (continuum) filters. The resulting Lya 
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images are the product of careful modeling of both the stellar and nebular continua, facilitated 
by supporting HST imaging at A w 2200, 3300, 4400, 5500, Ha and 8000A, combined with Star- 
burst 99 evolutionary synthesis models, and prescriptions for dust extinction on the continuum. 
In all, the resulting morphologies in Lya, Ha, and UV-continuum are qualitatively very different 
and we show that the bulk of Lya emerges in a diffuse component resulting from resonant scatter- 
ing events. Lya escape fractions, computed from integrated Ha luminosities and recombination 
theory, are found never to exceed 14%. Internal dust extinction is estimated in each pixel and 
used to correct Lya fluxes. However, the extinction corrections are far too small (by factors from 
2.6 to infinity) to reconcile the emerging global Lya luminosities with standard recombination 
predictions. Surprisingly, when comparing the global equivalent widths of Lya and Ha, the two 
quantities appear to be anti- correlated, which may be due to the evolution of mechanical feed- 
back from the starburst. This calls for caution in the interpretation of Lya observations in terms 
star formation rates. The images presented have a physical resolution 3 orders of magnitude 
better than attainable at high redshifts from the ground with current instrumentation and our 
images may therefore serve as useful templates for comparing with observations and modeling of 
primeval galaxy formation. We therefore provide the reduced Lya, Ha, and continuum images 
to the community. 

Subject headings: galaxies: ISM — galaxies: starburst — ultraviolet: galaxies — galaxies: individual 
(Haroll, SBS 0335-052, IRAS 08339+6517, Tol65, NGC6090, Tololo 1924-416) 



1. Introduction 

The search for galaxies experiencing their first 
phase of violent star formation is one of the 
great challenges for observational cosmology. In 
this context, the importance of the Lyman-alpha 
emission line (Lya) can hardly be overstated. 
Young galaxies with low dust content should pro- 
duce strong Lya emission with equivalent widths. 



W hva , larger than IOOA (Chariot fc Fall 



1993 



Schaerer|2003| . With a rest wavelength of 1216~A 
the Lya line is observable from the ground in the 
optical window at redshifts z ~ 2 to 6, and in 
principle in the near IR as far as z ~ 20. 

Significant samples of galaxies can currently be 
established up to z > 5 using continuum (Lyman 



break, e.g. |Iwata et al.||2003| Stanway et al.| [2003 
Bremer et al.|2004||Shimasaku et al.|2005| or emis- 



sion line (Lya, e.g. Rhoads & Malhotra 200 1| Hu 
et~aLl[20U4l |Venemans et al.||2004| |Ouchi et aT 



2005) techniques. However, while there is now 
much information regarding the continuum mor- 
phologies of galaxies and their evolution with red- 



shift (e.g. [Windhorst et al. 2002 L otz et al.||2006 
Ravindranath et al. 2006 Law et al.|2007[ ), the sit 
uation is very different for Lya where morpholo- 
gies are all but unknown. At high-z (> 2) when 
Lya is observable from the ground, even the best 
seeing conditions result in observations with a spa- 
tial sampling of several kpc. At z < 2, space-based 



observation becomes a requirement, although still 
very few instruments have been capable of emis- 
sion line imaging in the UV. In the current study 
we take a first step towards remedying the almost 
complete lack of information regarding Lya mor- 
phologies of local star forming galaxies. 

First theoretical models proposed that primeval 
galaxies should be strong UV emitters both in 



continuum and line emission ( 


Partridge & Peebles 


1967 


Meier 


1976 


Sunyaev et al. 


1978 


), although 



early observations proved disappointing. However, 
at these times no high-z normal galaxies were 
known and the UV domain was terra incognita. 
An early attempt to study Lya emission from 
three local star forming galaxies using the Interna- 
tional Ultraviolet Explorer (IUE) found that Lya 
emission was only present in the most metal-poor 
galaxy studied (TMeier & Terlevichl 119811). Sub- 



sequent studies of larger samples (Chariot & Fall 



1993) confirmed a weak correlation of Lya and the 
metallicity of the gas; not an unpredicted result 
given that the dust content generally correlates 
with metallicity. However, corrections for dust 
were found not to reconcile the Lya /H/3 line ratios 
with those expected from recombination theory 



(Giavalisco et al. 19961 and it became clear from 
these early studies that some other processes must 
be at work. Indeed, resonant scattering of Lya 



in Hi had long been studied theoretically (Oster 
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brock||1962| |Adams|[l972| |Harrington||1974[ ). The 
sub-recombination line ratios and low Wh ya can be 
explained by multiple resonant scatterings locally 
trapping Lya photons, redirecting them and caus- 
ing them to traverse much greater integrated path 
lengths than non-resonant continuum or emission 
lines. Thus, even minute quantities of dust may 
result in significant or complete attenuation of Lya 
QNeufeld|p90| . 

A major breakthrough in the detection of signif- 
icant numbers of LAEs at high- 2 Lya came with 



Cowie & Hu (19981, and today Lya emission is 
arguably the most successful method for examin- 
ing the galaxy population at the highest redshifts. 
Probing intrinsically fainter objects than the Ly- 
man break technique, Lya emission serves either 
as the spectral signature through which to identify 
high-z galaxy candidates or the definitive redshift 
identification of candidate objects. The impact on 
astrophysical cosmology has been enormous with 
Lya surveys being used to study star-formation 



rates (Kudritzki et al. 2000 Fujita et al. 2003 



Gronwall et al.|2007 |, large-scale structure (Vene 



mans et al.||2004| |Ouchi et al.|[2005] |Nilsson et~aT 



2007 ) , test the ionization fraction of the intergalac- 



tic medium (Kashikawa et al.[2006 Dijkstra et al 



and identify potential sites of population 



ill star formation ( 


Malhotra & Rhoads||2002 


Shi- 


masaku et al.||2006 


1. 



While it's difficult to argue against the im- 
pact of the high-z Lya studies, some inconsis- 
tencies call for caution. Most notably that star- 
formation rates (SFR) in individual galaxies are 
often found to be lower when estimated by Lya 



when compared with UV continuum (Ajiki et al. 



2003 Taniguchi et al. 20051. Underestimates of 



the cosmic SFR density also tend to be under- 
estimated by larger values when measured from 



LAEs (Rhoads et al. 2003 Kodaira et al. 20031. 



Moreover, LAEs have an inherent selection bias 
since only sources which show emission are in- 
cluded. In the sample of ~ 800 UV selected 
galaxies at z ~ 3, (also known as Lyman Break 
Galaxies (LBG) |Shapley et aT1|2003| ) found a me- 
dian W^ya of ~ OA, with only ~ 25% showing 
WY,y a > 20A. How these values evolve towards 
fainter continuum levels is unknown. Furthermore 



( |Ouchi et "aT]|2003[ ) noticed that bright LAEs ap- 
pear more strongly biased against the underly- 
ing dark matter halo distribution than continuum- 



bright sources. While the semi-analytical models 
for the space density of LAEs seem to be con- 



verging with the observed values (Le Delliou et al. 



2005), uncertainties arc still present and therefore 
open to discussion. In order to derive robust quan- 
tities from the LAE population, it is vital to un- 
derstand the origin of these biases. 

Again, results from low-z imply that the inter- 
pretation of Lya observation at high-z may not 
be so clear-cut. High spectral resolution studies 
with HST showed the Lya properties of local star- 
bursts to be complex, in realizing that a low dust 
content is by no means a guarantee for observing 
Lya in emission. Again, little correlation is found 
between Lya emission or absorption and dust or 



metallicity ( 


Kunth et al. 1994 


Lequeux et al. 1995 


Thuan & Izotov|1997 


)• 


Kunth et al. 


(1998 


) clearly 



demonstrated that when Lya is seen in emission 
it is systematically P Cygni in form and associ- 
ated with signatures of an outflowing neutral ISM, 
while starbursts showing Lya absorption features 
show no such outflow signatures. These results, 



also noticed at high-z in LBG spectra (e.g. Shap- 



ley et al. 20031, have profound implications for 
the escape of Lya from starbursts: they suggest 
an outflowing medium to be a requirement in or- 
der for photons to avoid resonant trapping in the 



neutral medium. Hydrodynamic models (Tenorio- 



Tagle et al. 1999 1 predict an evolutionary sequence 
for the Lya profile as a function of the mechani- 
cal energy return from the starburst, demonstrat- 
ing how absorption, emission and P Cygni phases 
may all arise during different phases. |Mas-Hesse| 
et al. (2003) used this physical evolutionary sce- 



nario to interpret the range of Lya profiles ob- 
served at low-z as an evolutionary sequence in the 
starburst population. To get a feeling for the as- 
trophysics involved, we encourage the reader to 



look at Fig. 18 in Mas-Hesse et al. (2003) 



Much theoretical effort has gone into the study 
of Lya line profiles. Ahn (20041 and Verhamme 



et al. ( |2006[ ) have examined the formation of 
P Cygni profiles from Lyman break-type galax- 
ies, finding secondary emission peaks redwards of 
the P Cygni emission feature resulting from one 
or more backscatterings in the expanding neutral 



medium. Ncufcld (1991) and Hansen & Oh (2006) 



have examined the transfer of Lya in a multi- 
phase ISM. For models in which dust is embedded 
in dense Hi clouds, they find that Lya photons 
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can in fact preferentially avoid dusty regions and 
escape having experienced lower total dust optical 
depths than non-resonant photons. Thus equiva- 
lent widths may be artificially boosted, potentially 
causing ordinary star-forming objects to appear as 
more exotic systems. 

What the low-z spectroscopic studies have not 
been able to accomplish is examination of global 
values relating to Lya: spectroscopic observa- 
tions miss the fraction (often the majority) of 
light which comes from outside the aperture or 
slit. This is highly likely to be the case for Lya 
where resonant scattering may cause decoupling 
of Lya from other wavelengths and Lya photons 
may scatter in the galaxy and emerge far from 
their sites of production. Ionized holes or porosity 
in the ISM may also increase the chance of Lya es- 
cape along certain paths. Ultimately there is little 
evidence for the spatial correlation between Lya 
emission and UV continuum at all. These con- 
siderations were the motivation for our pilot low-z 
imaging study of six local star- forming galaxies us- 
ing the Advanced Camera for Surveys (ACS) on- 
board HST. The intent was to map Lya morphol- 
ogy and study the effects of dust, kinematics, and 
luminosity on the escape of Lya and the sample 
was hand-picked to reflect a range of these param- 
eters. First results were presented in |Kunth et al.| 
(120031). 



Besides Lya, images have been acquired in Ha 
and six continuum bandpasses between 1500 A and 
the I-band, enabling not only a comparison be- 
tween Lya and Ha, but also detailed modeling of 
the stellar population and dust. Detailed studies 
have been presented for two galaxies in the sam- 
ple, both low-metallicity UV-bright starbursts, 



ESQ 338-04 and Haro 11 (Hayes et al. 2005 2007 1. 



These studies found that emission and absorption 
varied on small scales in the central starburst re- 
gions, exhibiting very little correlation with the 
morphology probed by stellar continuum observa- 
tions. Damped Lya absorption was found in front 
of many UV-bright, young star clusters that must 
be a significant source of ionizing photons. Using 
Ha, it was found that when Lya does emerge from 
central regions, it is typically at surface brightness 
well below the value predicted by recombination 
theory. In contrast, large low surface brightness 
Lya halos were found surrounding the central re- 
gions where typically Lya/Ha ratios and equiva- 



lent widths exceed those predicted for recombina- 
tion. Furthermore, modeling of the stellar popu- 
lation revealed that Lya appears to emerge from 
regions either too old to produce ionizing photons 
and/or too dusty to permit the transfer of Lya. 
Synthesis of these results points to the resonant 
scattering process being highly significant in the 
morphological structure. 

It should be noted that these HST imaging ob- 
servations are the only ones in existence that allow 
the comparison of Lya with Ha and stellar contin- 
uum and that they have a physical resolution, in 
the most distant case, 2 orders of magnitude bet- 
ter than that attainable at high-z under the best 
observing conditions. These data have now been 
fully processed, and continuum subtracted. De- 
tails of the continuum subtraction methodology 



can be found in Hayes et al. ( 2008 ) , hereafter re 



ferred to as Paper I. In this article we discuss the 
Lya line morphology, comparing it against Ha and 
stellar continuum, and compute integrated frac- 
tions of escaping Lya photons. The resolution of 
the data is so high that we believe them to be 
of interest for comparison with high-z samples of 
Lya-emitting galaxies and also for the computa- 
tional modeling of Lya radiative transfer. As a 
service to the community we release the processed 
data on the Centre de Donnees astronomiques de 
Strasbourg ( CDS). 

The article proceeds as follows: in Sections [2] 
and [3] we present the observations and data pro- 
cessing methods; in Section [4] we present the re- 
sults, which we discuss in Section [5j Section [6] is 
dedicated to future perspective and possibilities; 
in Section [7] we present our concluding remarks 
and describe how to retrieve the released data. 

Throughout this paper we will adopt Hq = 
72 km s^ 1 Mpc -1 for the Hubble constant. 

2. Observations 

Six galaxies were selected for observations with 
HST. Since this is a pilot study, our sample is 
hand-picked. We aimed to provide a reasonable 
coverage of the starburst luminosity, metallicity 
and dust content. Four known Lya emitters were 
included from the 
(19921 and 



ncy 



IRAS 08339+6517, NGC6090, 



samples of |Calzetti fc Kin 



Kunth et al. 



1998| ) : Haro 11, 
and ESQ 338-04; 



and two non-emitters from Thuan & Izotov ( 1997 1 
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SBS 0335-052 and Tol 65. Since Lya regulation is 
affected by ISM kinematics, selection of emitting 
targets was also based on the kinematic informa- 
tion imprinted on the observed spectroscopic line 

Our sample is 



2003) 



profile (Mas-Hesse et al. 
summarized in Table [T] 

The targets were observed with the Advanced 
Camera for Surveys (ACS) onboard the Hubble 
Space Telescope under two general observer pro- 
grams: GO 9470 for acquisition of far UV data 
in the F122M and F140LP filters using the solar 
blind channel (SBC); and GO 10575 for obtain- 
ing Ha, near-UV, and optical broad-band observa- 
tions with the High Resolution (HRC) and Wide 
Field (WFC) cameras. The F122M observations 
were obtained during the 'SHADOW part of the 
orbit in order to reduce the geocoronal Lya back- 
ground. In addition some complementary archival 
WFPC2 data was used. The observations are pre- 
sented in Table [2] 

3. Data reduction and processing 

AlQ pipeline calibrated data from the dif- 
ferent cameras were drizzled to the same scale 
(0.025"/pixel) and orientation using the task MUL- 
tidrizzle in pyraf. Remaining pixel shifts 
were fixed using the GEOMAP and GEOTRAN, and 
any remaining cosmic rays were removed using 
credit. Next, all the images for a particular 
target were convolved to a common point spread 
function (PSF) using the image with the broadest 
PSF as a reference. This in general was F550M, 
except for SBS 0335-052 and Tol 65 where the 
archival I-band data had been obtained with the 
WF3 chip of WFPC2. Each image was then sky 
subtracted using fit/flat _SKY in midas using a 
low order (0 or 1) polynomial and a common sky 
region for the different passbands. 

Obtaining Lya line-only images has been shown 
to be a rather involved process, requiring the use 
of spectral model fitting, and is presented in Sec- 



tion 3.1 Model parameters selected for the fitting 
are discussed in Section 3.2 Fitting in turn re- 



quires a threshold signal-to-noise ratio to yield 
tolerable results, requiring binning of the data, 
discussed in Section l3~3l 



3.1. Continuum subtraction 

Going from F122M and F140LP to continuum 
subtracted Lya images is non-trivial for a num- 
ber of reasons: the on-line filter is quite broad 
with a red wing, the effective wavelength of the 
off-line filter is rather far to the red side of the 
Lya line, neutral hydrogen in the Milky Way ab- 
sorbs part of the flux in the F122M filter, and the 
UV continuum is less well described by a power 
law near Lya than at slightly longer wavelengths. 
In effect, adopting a simple assumption of a UV 
power-law (f\ oc A' 3 ) and neglecting the other ef- 
fects mentioned can produce results ranging from 
strong Lya emission t o damped absorption s imply 



by varying (3. Indeed, Hayes & Ostlin 



(2006) have 



shown that also for high-z imaging observations, 
estimation of the continuum at Lya may be non 
trivial. 

To overcome these limitations we have de- 
veloped a method for continuum subtraction of 
ACS/SBC Lya images. It is described in detail 
in (Hayes et al. 2008 hereafter Paper i) where 



we also perform thorough tests of the method- 
ology and parameter dependencies. In short, the 
method utilizes Starburst 99 ( |Leitherer et al.| 1999 



^-With the exception of Tol 65 due to guide star problems, 
see Sect. 4.5 



hereafter SB99) spectral evolutionary synthesis 
models and a set of UV and optical continuum 
band passes plus Ha to perform a spectral energy 
distribution (SED) fit of up to two stellar and one 
nebular components, with reddening treated as 
a free parameter. The output is the Continuum 
Throughput Normalization ( CTN) factor which 
determines how to scale the F140LP flux in or- 
der to accurately continuum subtract the F122M 
images. 

The CTN is primarily dependent on the age and 
reddening of the stellar population dominating the 
UV flux. By sampling the UV spectral slope long- 
wards of 1500A and the 4000-break a non degener- 
ate solution for CTN is obtained. Additional op- 
tical bandpasses and Ha are included to constrain 
the contribution from an underlying older stellar 
population and ionized gas: Determining the con- 
tinuum level at Lya depends on fitting the stellar 
population age, which is sensitive to the strength 
of the 4000A break. Contamination of this feature 
by an aged population that does not contribute to 
the continuum at 1216 A could therefore result 
in erroneous age fits and estimates of CTN. Hence 
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Tabic 1: Target summary 



Target 


Alternative 


RA(2000) 


Dec(2000) 


E(B — V)mw 


log(nm)MW 


v r 


12+ 


M B 


Em/Abs a 


Ref 


name 


name 










(km/s) 


log(0/H) 








Haro 11 


ESO 350-38 


00:36:52.5 


-33:33:19 


0.049 


20.4 


6175 


7.9 


-20 


Em 


1 


SBS 0335-052 


SBS0335-052E 


03:37:44.0 


-05:02:40 


0.047 


20.6 


4043 


7.3 


-17 


Abs 


2 


IRAS 08339+6517 


PGC 024283 


08:38:23.2 


+65:07:15 


0.092 


20.6 


5730 


8.7 


21 


Em 


3 


Tol65 


ESO 380-27 


12:25:46.9 


-36:14:01 


0.074 


20.7 


2698 


7.6 


-15 


Abs 


4 


NGC6090 


Mrk 496 


16:11:40.7 


+52:27:24 


0.020 


20.2 


8785 


8.8 


21 


Em 


3 


ESO 338-04 


Tol 1924-416 


19:27:58.2 


-41:34:32 


0.087 


20.7 


2832 


7.9 


-19 


Em 


1 



a Em: known Lya-emitter, Abs: known Lycv-absorber 

Note. — E(B— V)mw i s the Galactic extinction according to Schlcgcl ct al. ( 1998| an d logfrnrrlmw is the Gala ctic Hi column 



density alon g th e line of sight. Reference refers to nebular metal licity measurements: 1: Bergvall . 
|et al.| ( |2006} , 3: |GonzaIez Delgado et aT] l |1998| ), 4: |Izotov et al.| | |2001| | 



Ostlin 



120021, 2 



Papadcros 



the need for a two population SED fit. The philos- 
ophy is similar regarding the B aimer jump which 
may be significant where Hn shells or filaments 
are well resolved, and is accounted for by treating 
nebular continuum independently 

The Ha images have been continuum sub- 
tracted within the SED fitting software and have 
been corrected for the contribution from [Nil] us- 
ing literature data. 

For all targets but two we have used two stel- 
lar plus one nebular component in the continuum 
subtraction (this is referred to as method v in Pa- 
per I) . For IRAS 08339+6517 and Tol 65 we did not 
have sufficient data for this purpose and used an 



law. For all galaxies we have adopted the SMC 



SED fit with fewer components (see Sect. 4.5 1 



3.2. 



Starburst 99 and continuum modeling 
parameters 



For this paper we use a standard setup of 



model parameters: A Salpeter (19551 IMF (a — 
2.35, dN oc M~ a dM) with lower and upper mass 
limits of 0.1 and 12OM , respectively. To model 
the starburst and underlying stellar populations 
we use instantaneous burst (single stellar popu- 
lations); Geneva tracks with enhanced mass loss, 
and Pauldrach/Hillier atmospheres for the UV 
were used throughout. While derived ages and 
reddenings are more sensitive to assumed model 
parameters, the derived CTN and hence Lya flux 



is indeed much less sensitive (see Hayes et al. 2005 
2008). In general, the most important concerns 



reddening law (Prevot ct al. 1984) since the intense 



UV radiation field of these galaxies tends to de- 
stroy the graphite component of the dust, leading 
to smooth extinction laws without the 2175 A fea- 
ture, and thus providing consistently improved fits 
over other laws (e.g. Mas-Hesse & Kunth 19991. 



The effect of changing this is discussed in detail in 
Paper I. 

For metallicity, we opt for the metallicity 
in Starburst 99 closest to the observed nebu- 
lar oxygen abundance (see Table [T|). Hence 
for IRAS 08339+6517 and NGC6090 we used 
Z = 0.02, for Haro 11 and ESO 338-04 we used 
Z = 0.004, and for SBS 0335-052 and Tol 65 we 
used Z — 0.001. Using Monte Carlo simulations 
(Paper i) we show that the accuracy of the contin- 
uum subtraction is not greatly affected, provided 
the assumed metallicity is within 50% of the true 
value. For all galaxies and metal abundances pre- 
sented here, that is indeed the case. 

3.3. Binning 

For the continuum subtraction to yield photo- 
metrically reliable results a signal-to-noise (S/N) 
in each passband of between 5 and 10 is required 
( Hayes et al.|2 008). This requires some smoothing 
or binning of the data in the low surface bright- 
ness regions. In order to simultaneously preserve 
the spatial information in the high S/N central re- 
gions, adaptive binning was used. We here make 



are the use of correct metallicity and extinction 



use of the Voronoi binning code of Diehl fc Statler 
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Table 2: Observations: HST filter names, central wavelengths, and exposure times in seconds. 



Target / Filter F122M F140LP F220W F330W F435W F550M FR656N F814W 
^central (A) 1274 1527 2257 3362 4312 5582 variable 8121 



Haro 11 


9095 


2700 


1513 


800 


680 


471 


680 


100 


SBS 0335-052 


9000 


2700 


1660 


800 


680 


430 


680 


4400" 


IRAS 08339+6517 


9000 


3000 


1590 


730 


360 c 


560 c 


1024 c 




Tol 65 


9095 


2700 


1728 


800 


680 


498 


680 d 


3400" 


NGC6090 


9095 


3000 


1683 


800 


680 


645 


680 


100 


ESQ 338-04 


9000 


3000 


isoo 6 


900 6 


800 b 


830 


1214 


520 6 



"Archival WFPC2/WF3/F791W image originally obtained in program 5408 (PI: Thuan) 

6 Archival WFPC2/PC data (F218W, F336W, F439W & F814W) from program 6708 (PI: Ostlin) 

c Data obtained with ACS/HRC 

d Data not useful 

e Archival WFPC2/WF3/F814W image originally obtained in program 6678 (PI: Thuan) 

NOTE.— Unless noted otherwise the F220W, F330W, F814W were obtained with ACS/HRC; and F435W, F550M, FR656N 
with ACS/WFC 



( 2006 ) which is a generalization of the binning al- 
gorithm of |Cappellari fc Copin ( 2003 1 . We use 



the F140LP image and its variance map to bin to- 
gether pixels to conglomerate groups (known as 
'spaxels') to improve S/N locally, while conserv- 
ing surface brightness. Each spaxel accretes pixels 
until a S/N of 10 has been attained or until the 
maximum spaxel size of 1600 pixels (1 □") was 
reached. The same binning pattern was then ap- 
plied to the other passbands. 

3.4. Uncertainties 

The accuracy of the continuum subtraction as 
a function of the input parameters is addressed in 



Hayes et al. ( 2008 1 , where noise is added to ideal- 



ized observations using Monte-Carlo simulations. 
In that article we consider only idealized, model 
data. These Monte Carlo simulations found the 
uncertainty (including the statistical uncertainty 
arising from noise in our fitting procedure, and 
systematic uncertainty accounting for our own es- 
timation of model parameters) on the Lya flux in 
a given spaxel to be a function of equivalent width. 
In general this has been shown to be below 30% for 
EW(Lya)=10 A and S/N=10 (the fiducial mini- 



mum value set by the binning, Section 3.3), with 



the error decreasing with increasing EW(Lya) and 
S/N. 

In this paper we have in addition performed 
Monte-Carlo simulations to assess the statistical 



uncertainties on the real data: to each spaxel, ran- 
dom errors have been added in accordance with 
the photometric uncertainties for each passband. 
For each galaxy 100 Monte Carlo realizations were 
performed. Thus, the 8 individual science im- 
ages for each object were regenerated 100 times 
with each spaxel randomly replaced by a Gaus- 
sian variate of the original, in accordance with 
its error. From these we created 100 new contin- 
uum subtracted images for each target, retaining 
all the important statistics of the spaxels in the fi- 
nal frame. We conclude that the statistical errors 
on the integrated Lya properties in Table [4] are 
very small, ranging from 4% in the case of Tol 65 
to 0.5% for ESO 338-04. The statistical S/N for 
each individual spaxel is of course smaller in gen- 
eral: the luminosity weighted averages range from 
(S/N) L = 0.26 for Tol 65 to 4.9 for Haro 11. Sta- 
tistical uncertainties on Ha properties are in gen- 
eral orders of magnitude smaller. 

A larger source of uncertainty is likely to be er- 
rors that result from sky subtraction and flat field- 
ing. The uncertainty of the sky level were, for each 
passband and galaxy, again fed into a Monte-Carlo 
simulation in a similar manner to that described 
above. A very conservative result is that flatfield 
and sky subtraction residuals can at most cause 
global errors on the 20% level. Since these kind 
of errors, when present, tend to have low spatial 
frequency they are likely to affect the global levels 
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rather than causing pixel to pixel variations. 

In addition we have estimated the uncertain- 
ties by calculating the standard deviation of the 
background of the resulting Lya level outside the 
galaxies, at a scale of 1 □", i.e. the maximum 
spaxel size (see Table [3]) . 

4. Results 

Various analyses are performed on the pro- 
cessed images: direct photometric analysis and 
comparison with model-determined properties. 
We examine both small-scale, local relationships 
and global quantities, and our results are pre- 
sented in a number of Figures and Tables. More 
detailed analysis for individual targets have ei- 
ther been presented previously ( Hayes et al.[2 005 



2007) or will be the scope of subsequent papers. 



4.1. Morphologies 

In Figs, [I] and [2] we show images of our targets 
in the F122M and F140LP filters, and continuum- 
subtracted Lya and Ha images. In this article 
we make use of adaptive smoothing (using fil- 
ter/ ADAPTiv in midas) and adaptive binning 



(Diehl & Statler 2006) techniques, both of which 



are widely used in the X-ray community. Smooth- 
ing results in images that are more pleasing to 
the eye, although less photometrically secure. We 
here show Lya images that have been obtained 
by continuum subtracting unbinned data, and 
an adaptively smoothed version of the same im- 
age. F140LP FUV contours are overlaid on the 
smoothed Lya image to guide the eye. Binning, 
on the other hand, strictly conserves local surface 
brightness and has consequently been used for the 
photometric analysis (Sect. 4.3 - 4.4 and Tables 
4 and 5). 

In Fig. [3] we show RGB color composites for 
our six targets, with Ha in red channel, F140LP in 
green, and adaptively smoothed Lya in blue. The 
intensity scale is set to show detail and in general 
the Ha and optical emission can be traced to much 
larger radii than here apparent. Both greyscale 
and false-color composite images show Lya mor- 
phologies that are strikingly different from those 
in UV and Ha. In all galaxies we see a mixture 
of both emission and absorption although Lya ap- 
pears not to emanate from regions of highest FUV 
surface brightness and instead seems to be most 



frequently seen from regions fainter in the contin- 
uum. 

4.2. Considerations for quantitative Lya 
photometry 

Absorption of the stellar continuum, both in 
the Milky Way and internally, present a couple of 
caveats against accurate photometric analysis of 
the continuum-subtracted Lya maps. A review of 
these caveats is thus warranted. 

Ultraviolet aperture spectroscopy of two of our 
targets (SBS 0335-052 and Tol65) have revealed 



damped Lya absorption (Thuan & Izotov 19971 



This is confirmed by the Lya continuum subtrac- 
tions for both these objects. However, all obser- 
vations in F122M are also subject to Galactic Hi 
absorption and in cases where the redshift is small, 
internal and Milky Way features may be blended. 
While our continuum subtraction software does 
correct for the Galactic absorption, internal ab- 
sorption features may be disguised if the features 
overlap. This is indeed the case for Tol65 and, 
to a much lesser extent, SBS 0335-052. In effect 
this may lead us to overestimate the emerging Lya 
flux. 

In several cases when Lya has been detected in 
high resolution spectroscopy it is found to have 
a P Cygni profile due to outflows in the neu- 



tral ISM that promote the escape of Lya (Kunth 
et al.|[l998| . When the profile is unresolved, as is 
the case in low-resolution spectroscopy or narrow- 
band imaging, the absorption on the blue side may 
cancel part of the emission. It should be noted, 
however, that scattering alone does not lead to a 
reduction in flux, but rather a modified morphol- 
ogy and/or redistribution in frequency. Indeed, 
Verhamme et al. (2006) note that P Cygni pro- 



files can arise as a natural consequence of radia- 
tion transfer in a moving medium, with no net loss 
of photons. 

Whether moving neutral media will give rise 
to net emission or absorption not only depends 
on the Hi density and outflow velocity but also 
on the relative contribution of the line and UV 
continuum, which of course are well resolved at 
low-z. If the continuum is strong, as in front of 
UV-bright star clusters, then any absorption could 
potentially be strong as well, increasing the like- 
lihood of measuring net absorption. Indeed, in 



S 




Fig. 1. — Greyscale images. From top to bottom: Reduced F122M and F140LP; continuum subtracted 
Lya resulting from unbinned images but with the ^45,1500 = 23 mag/D" isophotal mask generated from 
binned data overlaid; adaptively filtered Lya with F140LP isophotes overlaid; and continuum subtracted 
unbinned Ha. The isophotal mask from the binned data is used later for quantitative photometry. The FUV 
contours show 25 mag/D" as the faintest isophote, with each subsequent isophote representing an increase 
by 1 mag/D". All intensity scaling is logarithmic. North is up, East to the left. Images are square cutouts 
with the length of one side in arcsec given in parentheses. From left to right: Haro 11 (20") , SBS 0335-052 
(10") and IRAS 08339+6517 (20"). In some images residual features from the ACS/SBC detector gap are 
visible. These areas have not been included in the photometry. 
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Fig. 2. — Greyscale images. As in Fig.[T] From top to bottom: F122M, F140LP, Lya, adaptively filtered Lya 
with F140LP isophotes overlaid, and Ha. From left to right: (size in arcsec in parentheses): Tol65 (10") , 
NGC6090 (20") and ESQ 338-04 (20"). The intensity scaling is logarithmic. North is up, East to the left. 
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Fig. 3. — Color composite images of the whole sample. The blue channel shows continuum subtracted 
Lya, adaptively smoothed for the sake of presentation. Continuum subtracted Ha is shown in red. 1500A 
continuum is shown in green. The intensity scaling is logarithmic and arbitrary and has been adjusted to 
show interesting details. Each image has a size of 15"xl5", except for SBS 0335-052 (5" x 5") and Tol65 
(7.5" x 7.5"). For Tol65 the red channel shows a mixture of H/3, [Oin] and Ha (see text for details). East is 
left, North is up. Low surface brightness Lya emission can clearly be seen surrounding compact UV point 
sources, particularly in ESO 338-04, Haro 11, and the South West knot of NGC6090. 
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a few cases we find local minima in Lya surface 
brightness in front of UV-bright clusters that re- 
sults from just this effect. On the other hand, in 
faint regions there is not much continuum to ab- 
sorb. While the absorption feature itself may be 
undetectable in this regime, the blue wing of the 
emission profile may still be absorbed. When dust, 
of unknown distribution is added, it is impossible 
to know whether the 'absorbed' photons from the 
blue wing in any given pixel are (i) absorbed lo- 
cally, (ii) scattered and absorbed elsewhere, (iii) 
scattered and emitted elsewhere, or (iv) scattered 
and redistributed in frequency to the red wing and 
emitted locally. Note that these effects will also 
be significant in typical high-z observations where 
the processes governing Lya emission must be the 
same. What the true Lya flux really is in the 
case of P Cygni-blended absorption and emission 
is a matter of definition. Here we naturally define 
Lya as the sum of emission and absorption. 

4.3. Quantitative global photometry and 
properties 

We have extracted 'global' broad-band and 
emission line fluxes using a mask which was made 
by cutting the Voronoi-binned F140LP image at 
a surface brightness of level of /x as, 1500 = 23 
mag/D", corrected for Galactic foreground red- 
dening. The global quantities presented in Ta- 
bles [3] and [4] have been integrated over these 
apertures. Table [3] also presents integrated far 
UV (A = 1527A) luminosities and UV contin- 
uum slopes (/?), derived from the F140LP and 
F220W images. In deriving we used the pivot 
wavelengths, given in Table [3] but the effect lve 
wavelengths of a filter will in turn depend on the 
spectral shape, i.e. (5. Hence the (3 values pre- 
sented are only approximate, although since they 
are not actually used in the analysis, this is good 
enough. 

Table [4] presents the derived Lya and Ha quan- 
tities. Regarding Lya, we firstly quote the mea- 
sured continuum-subtracted flux, /Lya, which is 
the sum of all pixels (positive and negative) within 
the aperture. In addition we quote the sum of only 
the positive pixels, f^ ya , thereby rejecting all pix- 
els that represent absorption of the stellar contin- 
uum. Note that our ability to resolve emission and 
absorption on small scales depends on the distance 
to the target. The integrated Lya flux within our 



apertures is in five cases positive and one case neg- 
ative. The Lya/Ha ratios are an order of mag- 
nitude lower than the recombination value of 8.7 



(Case B, Brocklehurst 19711. This indicates sig- 
nificant attenuation of Lya due to dust, possibly 
aided by radiative transfer effects from Hi. The 
Lya equivalent widths (W^ya) range from -37 to 
+39 A, i.e. significantly smaller than for the the- 
oretical predictions for young starbursts. To first 
order equivalent widths should be unaffected by 
reddening if the line and continuum are equally 
attenuated. The Ha equivalent widths (Who) ap- 
pear anti-correlated with W\, ya if anything (see 
Fig [7]). For pure recombination one would ex- 
pect Wn a and Wh ya to be positively correlated, 
so clearly other effects than ordinary dust redden- 
ing must be at work. 

In Table[5]we present quantities directly related 
to the escaping fractions of Lya and dust extinc- 
tion, integrated over the whole region of the mask. 
Firstly we present the Lya escape fraction which 
has been calculated from: 



esc, Lya 



EU(8-7x^W) (1) 



where fu a ,c is the Ha flux corrected for the in- 
ternal reddening found in our SED fit. This red- 
dening is actually that of the stellar continuum 
and may hence be an underestimate for the neb- 
ular gas (e.g. Calzetti et al. 2000), which may 
in turn cause us to overestimate fesc.Lya- See 



Atek et al. (2008) for an analysis of the nebular 
gas extinction in these galaxies. We find that 
/esc,Lya varies between (SBS 0335-052; actu- 
ally absorption so the measured escape fraction is 
negative) and around 14% (IRAS 08339+6517and 
ESO 338-04). f esc Lya 1S a ^ s0 presented, where 
only positive Lya pixels are included in the sum- 
mation. Rejection of negative pixels results in 
little change in the Lya escape fraction for the 
galaxies that show few or weak absorption centers, 
(e.g. IRAS 08339+6517), significantly increases 
it in some cases (e.g. Haroll and Tol65), and 
changes its sign in the case of SBS 0335-052, for 
which the escaping fraction of detectable photons 
is around 2%. In 
above 20%. 

Since we have an estimate of the internal red- 
dening, we deem it meaningful to also use this 
value to de-redden the Lya fluxes. This makes 
it possible to compute dust-corrected Lya/Ha 



no cases does /+ c Lyct reach 
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Table 3: Integrated broad band fluxes and magnitudes. 



Target r eqA500 A = 1527 2255 3360 4320 5580 8115 log(L F 

( ) />. fx fx fx fx fx 

(kpc) tuab "iAB toab toab '"ab i«ab 



Haroll 


5.4 


5.21E-14 


2.16E-14 


1.09E-14 


8.65E-15 


5.12E-15 


3.17E-15 


10.3 


-2.25 




2.2 


14.9 


15.0 


14.9 


14.6 


14.6 


14.3 






SBS 0335-052 


3.1 


1.24E-14 


4.38E-15 


1.71E-15 


9.71E-16 


4.08E-16 


1.95E-16 a 


9.3 


-2.66 




0.83 


16.4 


16.7 


16.9 


16.9 


17.3 


17.4 






IRAS 08339+6517 


5.9 


9.04E-14 


3.83E-14 


2.37E-14 


2.50E-14 


1.70E-14 




10.4 


-2.21 




2.2 


14.3 


14.4 


14.0 


13.4 


13.3 








Tol65 


2.3 


4.26E-15 


1.91E-15 


7.46E-16 


4.65E-16 


1.95E-16 


9.30E-17 b 


8.5 


-2.05 




0.42 


17.6 


17.6 


17.8 


17.7 


18.1 


18.1 






NGC 6090 


4.0 


1.71E-14 


9.86E-15 


6.81E-15 


6.84E-15 


5.18E-15 


3.60E-15 


10.1 


-1.41 




2.3 


16.1 


15.8 


15.4 


14.8 


14.6 


14.2 






ESQ 338-04 


7.6 


1.11E-13 




2.16E-14 


1.81E-14 C 


1.25E-14 


5.76E-15 b 


9.9 


-2.07 




1.4 


14.1 




14.1 


13.7 


13.6 


13.6 







a Wavelength 7880 
^Wavelength 8012 
c Wavelength 4310 

Note. — Integrated broad band fluxes and magnitudes within the isophotal mask (/iisoo,AB = 23 mag/arcsec 2 ). The second 
column gives the equivalent radius (r eq = \J Arca/7r) within which the fluxes were integrated in units of arcseconds (first line) 
and kpc (second line). First line gives fx in cgs units (erg s — 1 cm -2 A -1 ), second line gives magnitudes in the AB-system. 
Lfuv = A X L\ in units of Lq^oI = 3.8 • 10 33 erg/s, and A = 1527A. The UV-continuum slope is derived from /1500 vs. /2200 
except for ESO 338-04 where we used /3360 rather than the available F218W image which is not as deep and quite affected by 
CR-residuals. All quantities have been corrected for Galactic extinction using the Schlegel et al. (1998) prescription. 



ratios and dust-corrected Lya escape fractions. 
Of course, if standard recombination theory was 
directly applicable, the corrected Lya /Ha ra- 
tio should be 8.7 and the dust-corrected Lya 
escape fractions should be unity. Instead, the 
extinction-corrected Lya fluxes are well below 
standard recombination predictions with factors 
ranging from 2.6 to infinity. In Table [5] we present 
ft = E /Lya,c/(8.7 x ^/hq,c), i.e. the ratio 
of the extinction-corrected Lya/Ha to the re- 
combination value and equivalently refers to the 
extinction-corrected 'escape fraction' 

Comparing f eS c,hya with 1Z reveals that (for 
galaxies with net emission) the latter is a fac- 
tor 2 to 5 larger than the former. ESO 338-04 
shows a Lya escape fraction of almost 40% after 
de-reddening, and IRAS 08339+6517 approaches 
30%. The ratio of 1Z to / esc ,Lya is not significantly 
different when regarding positive pixels only, and 
therefore not presented. The emerging picture 
is hence that simple reddening corrections can't 
come close to explaining the found escape frac- 
tions, meaning that in addition to simple dust 
reddening, other effects, i.e. resonant scattering 



and ISM kinematics, appear to dominate the Lya 
transport. 

Of course, these Lya escape fractions would be 
even lower if we have underestimated the total ion- 
izing continuum radiation, i.e. if the true redden- 
ing of Ha photons is larger than found by our SED 
fitting which could happen if there were some very 
deeply embedded sources that are completely ab- 
sorbed in the UV. In Atek et al. (20081 we use 



ground based Ha and H/3 images to estimate the 
dust extinction and consistently find even smaller 
escape fractions. 

4.4. Small scale photometry and proper- 
ties 

In Figures [4] and [5] we show the Lya surface 
brightness vs. Ha, UV continuum, and 0. The 
Voronoi binned resolution elements used to mea- 
sure surface brightness vary in size between indi- 
vidual pixels in the central regions and 1600 pixel 
(ID") spaxcls in the faintest regions. Here we 
do not mask out the fainter regions as described 
for the integrated photometric measurements. In- 
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Tabic 4: Measured emission line properties. 



Target name 


/Lyo 








/Ho 






Lya /Ha 


Haro 11 


9.30E-13 


15.6 


8.40E+41 


4.0E-15 


2.48E-12 


704 


2.24E+42 


0.374 


+ only 


4.32E-12 


25.5 


3.90E+42 










0. 625 


SBS 0335-052 


-6.77E-13 


-36.8 


-2.59E+41 


2.9E-15 


3.61E-13 


1434 


1.38E+41 


-1.76 


+ only 


2.97E-H 


1.62 


1.12E+40 










0.078 


IRAS 08339+6517 


3.05E-12 


39.2 


2.37E+42 


3.4E-15 


2.37E-12 


109 


1.84E+42 


1.29 


+ only 


3.52E-12 


40.7 


2.73E+42 










1.54 


Tol65 


8.48E-14 


17.5 


2.82E+40 


2.3E-15 


1.85E-13 


1300 


3.41E+40 


0.458 


+ only 


1.47E-13 


15.5 


4.89E+40 










0.794 


NGC 6090 


4.20E-13 


28.3 


7.80E+41 


1.7E-15 


1.27E-12 


315 


2.36E+42 


0.330 


+ only 


5.71E-13 


38.2 


1.06E+42 










O.449 


ESO 338-04 


3.68E-12 


28.1 


6.89E+41 


3.7E-15 


2.61E-12 


570 


4.89E+41 


1.41 


+ only 


4.99E-12 


38.0 


9.35E+41 










1.91 



Note. — Integrated galaxy properties within the far UV isophotal mask (/iisoo, AB = 23 mag °~Lya is the 'spaxel— 

to— spaxel' standard deviation of the Lya background outside the photometric aperture and hence measures the fluctuations 
on a one arcsecond scale. Note that this does not entirely reflect the error from flatfielding and sky subtraction which is 
conservatively found to be smaller than 20%. All quantities have been corrected for Galactic extinction using the Schlegel et 
al. (1998) prescription. Fluxes and luminosities are in cgs units, and equivalent widths in A. The Ha luminosity and equivalent 
width of Tol65 has been derived from ground based observations (see text). The second line for each galaxy gives the result 



Lya vs. Ha: A general feature is a sharp 
boundary at maximum Ha surface brightness, at 
which Lya/ Ha ranges from unity and down- 
wards (a), except for SBS 0335-052 where Lya/ 
Ha reaches the recombination value. The large 
range of values (extending to negative values, but 
which are not shown in these log-plots) is either 
the result of dust reddening or resonant scattering 
decoupling the Lya photons from Ha. Another 
general feature is a ceiling in Lya surface bright- 
ness at ~ 10~ 13 erg/s/cm2/arcsec 2 (b). While 
Haro 11 and ESO 338-04 present similar features, 
we sec, in addition, a population of Lya-bright 
spaxels with Lya/ Ha at, or above, the recombi- 
nation value (c). In all galaxies an upper envelope 
is seen near or slightly above the recombination 
ratio, but at fainter Ha levels a tail with super- 
recombination values develops - an indication of 
resonant scattering. 

Lya vs. UV continuum: A rather sharp 
boundary is seen at UV-bright levels (~ 10~ 14 
to 10 -15 erg/s/cm 2 /D"/A), reminiscent of a star- 
burst intensity limit (d). From the combined UV 



when including pixels with positive Lya flux only. 

stead, for inclusion of a spaxel, the criterion of hav- 
ing reached S/N > 3 in the F140LP image must 
have been fulfilled. This allows us to probe the 
fainter surface brightness regions while maintain- 
ing data reliability. Resolution elements in which 
the desired threshold S/N > 10 (See |Hayes et al 



2008) has been met are shown in black, whereas 
pixels with 3 < S/N < 10 are shown in blue. 
Since the pixel distribution plots use bins of var- 
ious sizes, they cannot be used to estimate how 
the values arc globally distributed and instead, 
histograms are presented above and to the side 
of the scatter plots in order to give some feeling 
for this. Moreover, since the surface brightnesses 
are presented in log-space, only pixels showing 
emission are included. Dashed lines on the Lya 
vs. Ha plots show the recombination line ratio of 
fhya = 8.7x / Hq and also the case f Lya = / Ha . On 
the Lya vs. FUV plot the lines represent constant 
equivalent width of 1, 10, and 100A. These scat- 
ter plots show a number of characteristic features 
which we now discuss. For illustration several fea- 
tures are labeled (a-f) in the text and in Figures [2] 
and|U 
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Tabic 5: Derived and reddening corrected properties. 



Target fesc,L V a{+) K SFR Ha log(L FIR ) SFR H „/SFR FIR E(B-V), log(L FUViC ) 



Haro 11 


0.037 (0.062) 


0.17 


17.2 


11.1 


1.3 


0.067 


10.6 


SBS 0335-052 


-0.203 (0.009) 


-0.28 


1.1 


9.3 


2 


0.044 


9.45 


IRAS 08339+6517 


0.140 (0.168) 


0.28 


14.1 


11.0 


0.9 


0.034 


10.6 


Tol65 


0.053 (0.091) 


0.16 


0.26 


8.3 


5 


0.11 


8.94 


NGC6090 


0.027 (0.037) 


0.13 


18.2 


11.4 


0.4 


0.19 


10.8 


ESQ 338-04 


0.139 (0.189) 


0.39 


3.9 


9.8 


3.4 


0.076 


10.2 



Note. — f esc, Lya is the Lya escape fraction obtained by dividing the observed (Galactic extinction-corrected only) Lya flux 
by the Ha flux corrected for Galactic and internal reddening (as a measure of the total intrinsic Lya production) and the 
theoretical line ratio. TZ is the same quantity but using the Lya flux corrected for internal extinction. SFR^f a is based on the 
observed Ha flux in Table [4] (not corrected for internal reddening) and uses the calibration by Kennicutt (1998). log(Lpjn) 
is based on the 1.7 times the luminosity at 60/^m (Chapman et al. 2000) as detected by IRAS and listed by NED, except for 
SBS 0335-052 and Tol 65 which were not detected by IRAS and where we have instead used the Spitzer values reported by 
Engelbracht et al. (2008). The FIR SFR calibration in SFR// a /SFR P i R uses the calibration by Kennicutt (1998). E(B - V)* 
is the color excess found in the SED fit for the integrated fluxes within the same far UV isophote as above. 



and Lya peak towards fainter UV surface bright- 
ness levels we see an upper envelope, which starts 
at W-L,y a ~ 30A at the bright end and approaches 
1000A at fainter (~ 10- 17 erg/s/cm 2 /n"/A) lev- 
els where it forms a tail with W^ ya > 1000A (e). 
Again Haro 11 and ESO 338-04 present extra high 
surface brightness features with average equivalent 
widths of about 30 A. 

Hence, we see regions where line ratios and 
equivalent widths are significantly enhanced com- 
pared to what would be expected from recombina- 
tion theory and standard spectral synthesis mod- 
els. This may naturally be expected to some ex- 
tent for Lya vs. FUV since ionized regions may 
be spatially separated from the UV, but for Lya 
vs. Ha it must imply that some process, resonant 
scattering, is causing Lya to be emitted from re- 
gions where Ha is very faint indeed. 

Lya/ Ha vs. (3: Pixcl-to-pixel variations be- 
tween the Lya/Ha ratio and (3 again show little 
correlation as if Lya emission was independent of 
the local stellar reddening. A tail in the distribu- 
tion is also seen in these figures up to (3 > (f). 
These red UV slopes indicate that either the stellar 
population is too old to produce nebular emission, 
or the dust reddening is too great to permit the 
transmission of Lya. Either way, Lya is seen in 
emission from these regions, and a significant frac- 
tion of the pixels show ratios greater than that of 



case B recombination. 

UV surface brightness: In Fig. [6] we show 
the fraction of Lya fluxes emitted (or absorbed) 
within the masked regions as a function of the UV 
surface brightness. This is done by summing the 
flux in the pixels of the line-only image that fall 
within the surface brightness bins in the FUV im- 
age, then normalizing by the integrated flux within 
the masked region. Thus we can quantitatively as- 
sess how Lya emission/ absorption correlate with 
UV surface brightness. The same is also presented 
for Ha and UV continuum itself. In most cases 
Ha is clearly seen to be offset from the UV con- 
tinuum since, to some varying degree, the bright 
stellar sources and ionized nebulae can be resolved. 
While the bulk of the emitted Lya flux also ap- 
pears to emerge from regions of lower FUV surface 
brightness, the distribution is further offset from 
that of Ha. 

4.5. Comments on individual targets 
Haro 11 This target is extensively discussed in 



Hayes et al. (2007), where we also discuss X-ray 



data, the escape of Lya as a function of the age of 
the stellar population (field stars and stellar clus- 
ters) and nebular reddening (from ground based 
observations). |Kunth et al. ( 1998[ ) found an asym- 
metric Lya line in emission with no sign of Lya 
emission at restframe 1216 A, suggesting Lya es- 
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Fig. 4. — Pixel scatter distribution showing all 'spaxels' for which S/N > 3 in F140LP. Each spaxel, which 
may contain up to 1600 pixels of size 0.025" is represented by a single point. Black points show spaxels 
where S/N > 10 was attained, blue points show pixels in the range 3 < S/N < 10. For the histograms, 
however, each spaxel is represented by the number of pixels that contributed to it, in order to show the true 
distribution. Since the ordinate shows the logarithm of the Lya surface brightness, only spaxels with positive 
Lya flux are included. Error-bars show the median statistical error in each magnitude bin for the ordinates 
and abscissae (placed diagonally to avoid overlap). Left: Lya vs. UV surface brightness, here represented 
by the modeled continuum at 1216A, i.e. /iAi527 x CTN. The dashed lines show constant Lya equivalent 
widths of 1, 10, and 100 A. Center: Lya vs. Ha surface brightness, both in units of erg/s/cm 2 /arcsec 2 . The 
dashed lines show Lya = Ha, and Lya = 8.7Ha (Case B recombination). Right: Lya/Ha line ratio vs. UV 
continuum slope /3 as derived from the fluxes in F140LP and F220M. Top panels: Haroll. Middle panels: 
SBS 0335-052. Bottom panels: IRAS 08339+6517. 



16 



■ — 

f2 



- 
c 

y 
o 
z 



C 




lug (Lyra con*, surf .brigtrt . } 



-3 











1 II 1 1 A , 


■ f ^ 

' y * * 

'"-H * * * * 



] 



lug (Lyra £:[)nt. surf. bright.) lug (Hn tmrf.bright.) lug { N ) 



§ -13 



■It 










log ( N ) 



E 



w 1 




H 2 ■> I J 

j I rry ) 



: N ■ 



T--1 




-»'""'" ' f X 
^^^^^^^ 





log (Lycy colli. surf.brigtiL.} log (Ha surf, briglit.) log ( N ) 



3 2 . :i I ■= 
/ ( 1T\ \: V I 



log ( N ) 



Fig. 5.— Same as Fig. [3] but for Tol65 (top), NGC6090 (middle) and ESO 338-04 (bottom). For Tol65 we 
do not show Ha since we only have reliable data at ground based resolution. For ESO 338-04 /3 was derived 
from the F140LP and F336W data. See caption of Fig. [4] for further details. 



17 



(2.2 



0.1 



0.2 



0.1 



'-- 

u 
::. 
■■ 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- Harol 1 ■ 

✓ 

#*•*«£*:*. V 
..* ,* **• . . 
** * i * • - 


1 i i i 1 i i i 1 i i i 1 i i i 1 i ii | ii 
- E50 338 

* 

# \ A - 

***** 

-7 '- 

: a ■.* 

« 4k 
P ■ *' ■ * # 

•j* ». '*'*** - 


■ ■■■>'* 

~ i i 1 i i i 1 i i Jfcl 1 i i i 1 i i i 1 i i " 


' . — 

■V / 

~ 1 i i i 1 i i i 1 i i / 1 i i i 1 i ■ ■ 1 ■ ■ 


1 1 | 1 1 1 | 1 1 1 | 1 1 1 | 1 1 1 | 1 1 

NGC S090 

* v ■ 

; - .S *** **- 

* .ft" ■ - 


| i i i | i i i | i i i | i i i | i i i | i i 

■ IRAS 06+65 

■ ■ 

* J •mm 

■b 

■ ■■ 

: f 


t - 

~ i i 1 i i i 1 r i i 1 i i i 1 i i i 1 i i _ 


" 1 i i i 1 i i i; 1 i i i 1 i i i 1 i i i 1 i i" 


i i | i i i | i i i | i i i | i i i | i i _ 

SBSG335 
_ z 

ft**, 

: i» * * _ . 


_ | I I I | I I I | I I I | I I I | I l y_ 

7 Tol 65 ■ 

- 

■ 

■ 

: i .■ : 


,■■■»«■■■ 

■. * 

— i i 1 i i i 1 i i i 1 i i i 1 i i i 1 i r~ 


■ 

1 i i i 1 i i i l i i i 1 i i i 1 i i i l i — 



0.3 
0.2 
0.1 


-0.1 

-0,2 



1+ 16 18 SO 22 

Mfuv 



12 14 16 18 30 22 
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the /Mb, 1500 = 23 mag/D" mask) flux as a function of /iAB,i500 (corrected for Galactic reddening only) 
summed in 0.5 magnitude bins. Stars (green) show 1500A, filled squares (blue) show Lya, and asterisks 
(red) show Ha. Triangles (magenta) show the cumulative Lya escape fraction, determined from Ha - they 
always converge at the escape fraction presented in Table [5] Since the aperture mask is defined by making a 
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cape due to an outflowing neutral medium. In 
addition blueshifted metallic absorption lines in- 
dicated the presence of outflows along the line of 
sight. Bergvall et al. ( |2006 1 claimed a detection of 
Lyman continuum (LyC) photons escaping from 
this galaxy (but see also Grimes et al. 2007) and, 
since LyC does not scatter, porous ionized holes in 
the ISM are likely to exist along some sight-lines. 
The galaxy is very luminous in Ha and has a high 
inferred star formation rate, still no Hi or CO has 
been detected ( Bergvall et al. 2003[), the latter may 



be related to the low metallicity Bergvall fc Ostlin 



(2002) 



This galaxy consists of three major star form- 
ing condensations although Lya is seen directly 
in emission from only one (the NW knot; C in 
the nomenclature of |Kunth et al. | |2003[ ). It is 
this knot that causes the high surface brightness 
features discussed in Section |4.4| The remaining 
two knots (one of which dominates the Ha out- 
put) show Lya absorption centered on the clus- 
ters themselves, with a patch of faint diffuse emis- 
sion leaking from the condensation to the SW. The 
overall Lya flux escaping is dominated by a dif- 
fuse component with Lya/Ha ratios larger than 
the recombination value, indicating that most Lya 
leaves the galaxy after multiple scatterings (see 



Hayes et al. ( 2007 1 for a description of the Lya sur- 



face brightness profile). In the Lya-bright eastern 
knot Lya/Ha is close to the recombination value, 
and the EW is high, indicating that Lya emission 
from this region may be avoiding resonant trap- 
ping and escaping more or less unhindered. 

SBS 0335—052 Previous attempts have failed to 
detect anything but damped Lya absorption from 
this galaxy ( Thuan fc Izotov|1997 ). Here we show, 
that while the central region is indeed character- 
ized by Lya absorption there are regions from 
which Lya photons are emitted. Compared to the 
other galaxies, the fraction of Lya emitting spaxels 
is small, but significant at the 3cr level. Integrated 
over the galaxy, it is a net absorber of Lya. 

The brightest Lya emitting spaxels have Lya/Ha 
somewhat higher than the recombination value, 
which could result from resonant scattering, ISM 
geometry, or the production of additional Lya 
photons by collisional excitation. This galaxy ex- 
hibits a low surface brightness tail just as the 
brighter Lya emitters do. 



The lack of a net positive escape at any surface 
brightness level is somewhat puzzling as the Lya 
photons ought to get out somewhere, unless the 
density of static Hi is so large that even moderate 
amounts of dust will destroy Lya. 

IRAS 08339+6517 This is a spiral with a nu- 



clear starburst. Kunth et al. ( 1998) found P-Cygni 



like emission, and Mas-Hesse et al. (2003) found 



further support for an ISM-kinematic regulation 
of the Lya escape and evidence for recombinations 
in the outflowing shell. Lya/Ha again supports a 
mixture of dust and kinematic effects in the ISM. 

Since no F814W image was obtained for 
IRAS 08339+6517, our data-points are not suf- 
ficiently spread along the SED to fit two stellar 
populations, and continuum subtraction of Ha 
also becomes less reliable. For IRAS 08339+6517 
we therefore resort to a single component spectral 
fit using the standard stellar+nebular template 
from Starburst 99, i.e. method in in the terminol- 
ogy of Paper I. 

Tol 65 This is the faintest galaxy in the sam- 
ple, both in luminosity and apparent magnitude. 
In our HST Ha image, the galaxy happened to 
fall outside the monochromatic patch of the ramp 
filter and could not be used. Hence, we had no 
information about the distribution of ionized gas 
at HST resolution and could not model the nebu- 
lar emission as a separate component. Instead we 
used two stellar components where we used the 
option of Starburst 99 to include the nebular con- 
tinuum in the model (Leitherer et al. 19991. In 



addition, the ACS/SBC observations of Tol 65 suf- 
fered a guide star acquisition failure meaning that 
only one of the two guide stars were acquired, with 
the consequence that the telescope pointing could 
be drifting about the axis defined by the direc- 
tion to the single guide star used. The images 
were registered manually rather than using mul- 
tidrizzle since the astrometric keywords could 
not be trusted. We have carefully examined the 
SBC images and found no evidence for a signifi- 
cant telescope drift during the observations, but 
the results for Tol 65 should be viewed with some 
caution. 

To get a rough idea about the ionized gas dis- 
tribution we retrieved a WFPC2 F555W image 
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of exposure time 1600s from the archive. The 
F555W and F550M filters have very similar central 
wavelength but, whereas F550M does not transmit 
any strong emission lines, F555W transmits H/3, 
[Om]AA4959,5007, and partly Ha. Thus, by sub- 
tracting F555W-F550M, scaled by their PHOT- 
FLAM values, we get a map of the strong emis- 
sion lines, that can be used for qualitative pur- 
poses, e.g. in making the color image in Fig|3] 
For quantitative estimates of the Ha flux and 
equivalent width, we have used ground-based data 
obtained from ESO's New Technology Telescope 



(NTT) ( Atek et al.||2008|). Pr evious HST/GHRS 
spectra ( Thuan & Izotov 1997 ) have indicated that 



the brightest knot is a damped Lya absorber. The 
UV morphology of Tol65 is also rather complex 
but, assuming that the NW knot was acquired by 
the GHRS PEAK-UP maneuver, photometry per- 
formed in a 2" x 2" aperture consistently finds this 
to be an absorber. 

Overall we find that the continuum-bright re- 
gions arc dominated by absorption but we also see 
Lya in emission from fainter regions. Tol65 is 
the galaxy for which the data quality is worst and 
the uncertainties are largest. Due to the effects 
discussed in Sect. |4.2[ the positive net Lya flux 
observed from this galaxy may not be real. 

NGC 6090 This luminous starburst has a clear 
merger signature and is also very luminous in the 
IR. The star formation is concentrated to two re- 
gions both which show clear Lya emission with 
emission/ absorption varying on small scales. The 
Lya/Ha ratio scatters between above the recom- 
bination value and absorption. The common tail 
of Lya/Ha ratios is seen, extending to values that 
exceed the recombination value, and representing 
diffuse regions outside the star forming complexes. 

ESO 338-04 The Lya morphology of this galaxy 



was discussed in Hayes et al. (2005) using a less 



sophisticated method (similar to method in in the 
terminology of Paper I) for continuum subtraction 
which was based on a single stellar population 
and did not include a separate treatment of the 
nebular continuum. Our new results are, how- 



ever, in good agreement with Hayes et al. (20051. 



The starburst is dominated by a rich population 



of young massive star clusters ( Ostlin et al 



1998 



rich and detailed structure consisting of shells and 
filaments. In contrast the Lya emission is much 
more extended and clearly asymmetric with the 
main emission in the N-S direction, where there 
are signs of outflows from other investigations (e.g. 



Ostlin et al. 2001 



2007 ) . Central patches of Lya 



absorption are seen across much of the central star 
forming complex with a noticeable streak in the 
E-W direction. This is not seen in age or E(B—V) 
maps and is most likely the result of a filamen- 
tary Hi structure lying in front of the starburst. 
Unique for this galaxy is the clear vertical tail to 
super recombination values in the Lyavs Ha dis- 
tribution occurring also at bright levels (Knot A 
in the nomenclature of Hayes et al. (2005). See 
Fig. §. 

5. Discussion 

All targets reveal Lya morphologies qualita- 
tively very different from those revealed by the 
FUV continuum or Ha. Two Lya non-emitters 
(SBS 0335-052 and Tol 65) were on purpose in- 
cluded in the sample, but both show signatures 
of Lya leakage, and one of them is actually found 
to be a net Lya emitter (Tol 65), although this 
is our object with the poorest data quality. In 
the remaining cases (Haroll, IRAS 08339+6517, 
NGC 6090, ESO 338-04) where the Lya signature 
is very clear and well-structured, the equivalent 
widths are found to be low with line-ratios signif- 
icantly lower than the recombination value, and 
escape fractions well below unity. 

Our targets span luminosities from the dwarf 
regime to that of z ~ 0.15 compact UV- luminous 
galaxies (UVLGs) that have been proposed as lo- 
cal analogs of LBGs by virtue of their luminosity 
and morphology ( |Overzier et al. 



2008). Indeed, 



Haro 1 1 was used in their paper as a local com- 
parison object. Although their selection did not 
favour inclusion of galaxies with multiple nuclei, 
such as Haro 11 or NGC 6090 our four most lumi- 
nous galaxies otherwise have very similar proper- 
ties to UVLGs. Morphologically, and in terms of 
rest rest-frame colors and luminosities our sample 
galaxies are also similar to star bursts at z ~ 1 



studied in deep HST U-band images (de Mello 
etld]|2006l ). 



2003), around which strong Ha is seen following a 
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5.1. Emission scales and porosity 

In front of the star forming regions we find 
strong variations in both the Lya images and 
equivalent width maps, ranging between absorp- 
tion and emission. This phenomenon is seen on 
scales ranging between the smallest we can probe 
reliably (0.1", corresponding to a physical sizes 
ranging between 20 and 60 pc) up to sizes an or- 
der of magnitude larger. The scatter plots (Figs. 
[4] and |5| clearly show the large variation of Lya 
equivalent width (for positive spaxels) from the 
regions of peak UV and Ha surface brightness. 
There is no apparent characteristic size over which 
these variations are seen. That is, the Lya flux or 
Whya may be radically different, e.g. switching 
from absorption to emission, even on the smallest 
scales we can probe. This supports the model in 
which Lya may escape through favorable paths in 
a porous and inhomogcncous ISM. However, it is 
also clear that this is not the main mechanism by 
which Lya photons escape galaxies, because the 
net emission, if any, is always dominated by dif- 
fuse emission from regions with low UV and Lya 
surface brightness. 

In Figure [6] we have presented the fractional 
contribution to the total Lya and Ha emission as 
a function of UV surface brightness. Four out of 
five targets (Tol 65 was not observed in Ha) show 
a systematic offset between Ha and FUV (red 
lines appear shifted towards fainter levels com- 
pared to green), which can be expected for two 
reasons: Firstly, the ionization boundary will be 
extended away from clusters, where UV surface 
brightness is lower and may be resolved in the im- 
ages. Secondly, stellar mechanical feedback may 
blow out the ionized ISM in shells and filaments, 
further decoupling the nebular and stellar compo- 
nents. Likely this is true also for Tol 65 since, like 
in the other four cases, Lya is offset with respect 
to FUV. Lya is in general further (compared to 
Ha) offset towards lower FUV surface brightness 
levels, demonstrating the importance of additional 
scatterings affecting Lya. 

The only object that does not show a system- 
atic offset between FUV and Ha surface bright- 
ness is SBS 0335-052. |Mas-Hesse et all (p03| 



suggested that damped Lya absorption profile ob- 
served in this galaxy to be due to the youth of 
the starburst episode (Wn a = 1500A); insufficient 



time has elapsed to generate a strong wind and 
accelerate the surrounding Hi. The fact that we 
observe Ha and FUV emission from SBS 0335-052 
to be rather tightly correlated could also be a re- 
flection of this fact, and is consistent with the in- 
terpretation. The Lya line, on the other hand, 
almost completely mirrors that of the FUV, and 
is seen only in strong absorption in the central re- 
gions. This is naturally an effect of absorption of 
the stellar continuum as discussed in Section |4~2| 
where more flux appears to be absorbed at higher 
surface brightness. 

IRAS 08339+6517 and, to a lesser degree, 
NGC 6090 shows some degree of correlation be- 
tween the relative Ha emission and Lya. However, 
at intermediate UV surface brightness ((J,ab, 1500 = 
18 to 20 mag/D") systematically less Lya is 
emerging with respect to Ha than at fainter con- 
tinuum levels. Again this is indicative of the reso- 
nant scattering of Lya increasing typical physical 
scales of the emission and resulting in low-level 
emission outside the starburst region. 

The remaining strongly emitting galaxies, 
Haroll, and ESO 338-04 show strong and sys- 
tematic offsets between the surface brightnesses 
where the bulk of the Ha and Lya emerge. This 
is nicely illustrated also in Figures Tj2 and [3] 
where large scale Lya emission is seen from out- 
side the star forming regions. In Haro 11, the peak 
in the output Lya luminosity distribution is offset 
from Ha by about 2.5 magnitudes in [i ab,isoo- I n 
ESO 338-04 the bulk of the Lya output is shown 
clearly to be at faint UV flux levels although no 
clear maximum is seen and much Lya emerges 
from the faintest UV bin. 

Clear signatures of diffuse Lya emission, above 
the level of the background residuals is seen in 
both Haro 1 1 and ESO 338-04 at radial distances 
from the central star-forming region of around 3 
and 1.5 kpc, respectively. 

5.2. Lya escape fraction and dust 

Unlike non-resonant radiation, photon scatter- 
ing causes Lya to be emitted over larger areas 
than radiation in the adjacent continuum and at 
lower local signal-to-noise where it is more sus- 
ceptible to contamination by residuals from the 
sky-background and flatficlding. Thus the mea- 
sured Lya escape fraction may be a strong func- 
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tion of aperture which is evident from Fig [6j Had 
our data been deeper and we could have been 
more confident about photometry at very faint lev- 
els, the escape fractions may well have increased. 
For a case like SBS 0335-052, which is not very 
dusty, one would naively expect most of the Lya 
to emerge somewhere - indeed a significant frac- 
tion of Ha emission comes from outside our aper- 
ture mask. On the other hand, when we fit a Sersic 
profile to the Lya scattering component of Haro 1 1 



at good S/N (Hayes et al. 20071 and integrated 



the flux to an infinite radius, we obtained a flux 
consistent with aperture photometry within a few 
per-cent. For reference, to reach S/N=10 in a ID" 
aperture, a source of the same surface brightness 
as the average spaxel-to-spaxel sky variation (see 
Table Q would require over 30 orbits in SHADOW 
mode with SBC/F122M per target. 

Five out of six cases are found to be net Lya 
emitters with Lya escape fractions in the range 3 
to 14 percent. From the sixth galaxy (SBS 0335- 
052), we do detect some Lya emission leaking out- 
side the central starburst region. However, also 
at faint levels is the net integrated Lya flux nega- 
tive. It is possible that more Lya would have been 
detected in the case that our observations would 
have been deeper, allowing us to probe greater 
radii. This of course applies to all our galaxies: 
the regions for which the Lya and Ha fluxes have 
been integrated is significantly smaller than the 
total optical extent of the galaxies and we may 
be missing a non-negligible fraction of the diffuse 
emission. The Lya escape fractions should in this 
respect be seen as lower limits. On the other hand, 
some galaxies may contain star forming regions 
that are so dusty that not even any Ha escapes in 
which case the actual Lya escape fractions would 
be lower. 

Our modeling software provides us with an es- 
timate of the reddening of the stellar continuum 
below Lya; this is what is required for the subtrac- 
tion of the continuum from the F122M filter. How- 
ever, due to geometrical differences in the configu- 
ration of stars and nebular gas, stellar and nebular 



E(B—V) measurements are shown to differ i 


Maiz- 


Apellaniz et al. 


1998; 


Mas-Hesse & Kunth 


1999 


Chariot & Fall 


2000 


Calzetti et al. 


2000 


). Thus 



expected to minimize with increasing spatial res- 
olution of the observation, geometrical differences 
between the stellar and nebular components along 
the line of sight will always be present. Firstly, 
we know that Lya and B aimer lines must be pro- 
duced in the same ionized nebulae. Secondly, we 
may suspect that due to resonant scattering in- 
creasing the traversed path-length, Lya photons 
should experience at least the reddening of the 
Balmer lines, although certain geometrical config- 
urations may be envisaged that make this state- 
ment uncertain. Evidently, some circumspection 
must be exercised when comparing derived red- 
dening with the emission (or absorption) of Lya. 
With this cautionary remark in mind, we deem it 
interesting to de-redden the Lya fluxes (using the 
reddening found from the SED fit), recompute the 
Lya escape fractions, and explore to what extent 
pure dust attenuation can be said to regulate the 
fraction of escaping Lya photons. In that case we 
would expect to find 1Z = 1. After de-reddening 
we find Tl values between 13% (NGC6090) and 
40% (ESO 338-04). Hence, even in our dustiest 
case (NGC 6090) the observed stellar reddening is 
insufficient in explaining the attenuation of Lya. 

Figure [7] shows where the galaxies fall when 
f esc, Lya and 1Z are compared with E(B — V). 
Any correlation present for such a small sample 
must be viewed upon with skepticism, but the 
emitters with lowest E(B — V) do have higher es- 
cape fractions. ESO 338-04 is a dust-poor galaxy 
with the integrated measurement from this study 
finding E(B - V)= 0.076. Nevertheless, the 
applied dust correction does play a significant 
role, with the same conclusion being reached for 
IRAS 08339+6517, and also perhaps for Haro 11. 
1Z in some sense can be thought of as a measure 
of the importance of resonant scattering but the 
effect is so non-linear that quantification is not 
possible. 

To strengthen the analysis performed here we 
would need to independently map the nebular ex- 
tinction, at least mitigating the potential difficul- 
ties that arise from blended stars and gas. Such 
a study is performed at ~ 1 " resolution in (Atek 
et al.|2008 1 and we are currently exploring the use 



we have little information about the reddening in 
the gas in which Lya is formed. While differences 
in the respective reddening measurements can be 



of adaptive optics imaging in the Br7 line for two 
galaxies. 
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Fig. 7.— Scatter plots for the complete sample: Haro 11 (circle); SBS 0335-052 (triangle); IRAS 08339+6517 
(square); Tol 65 (diamond); NGC6090 (pentagon); ESO 338-04 (hexagon). Left: Lya escape fraction red- 
dening vs. E(B — V). Black points show the observed escape fractions, and red the escape fractions corrected 
for internal extinction. Since SBS 0335-052 is a net absorber, it is not shown. Right: Wi,y a vs. Wn a . Again, 
black points show the observed values, while for red symbols Lya and Ha have been corrected for the fitted 
E(B-V). 
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5.3. Equivalent widths 

Comparing Wt, ya or the Lya escape fractions 
with most of the parameters listed in Table [5] 
yields almost nothing in the way of trends. In- 
deed this study has never aimed to provide any 
statistical answers and with only a sample of six 
objects could hardly be considered statistically 
meaningful. Moreover, since the sample was de- 
signed to sample a range of starburst parameters 
and include emitters and non-emitters of Lya, any 
trends may purely be a feature of the selection. 
That said, the anti-correlation of W^ya and Who 
mentioned in Section |4.3| is indeed interesting as 
it seems to be rather tight, although this trend 
is not seen in the Lya escape fraction, with or 
without dust correction. From the Starburst 99 
models, we estimate that our lowest globally mea- 
sured Wft a corresponds to an age of around 6 Myr 
for an instantaneous burst, solar metallicity and 
Salpeter IMF. Higher Wn a naturally corresponds 
to younger ages, with equivalent widths of several 
thousand A expected during the first 3 Myr or 
so. In the studies of Tenorio-Tagle et al. (19991 



mote Lya escape. 



Oey et al. (20071 studied the ISM in nearby 



and Mas-Hesse et al. (20031 it is concluded that 



at the youngest times, Lya may be expected in 
absorption, followed by an emission phase when 
the mechanical energy return from star formation 
has accelerated the ambient medium. How the ex- 
pected equivalent widths of Lya and Ha co-evolve 
is difficult to predict, and instantaneous burst as- 
sumption is certainly an oversimplification when 
dealing with global properties. Nevertheless, what 
we are seeing may be a direct consequence of evo- 
lution. 

Related to the Lya escape fraction and equiv- 
alent width is the escape fraction of Lyman con- 
tinuum (i.e. ionizing) photons. Some Lyman con- 
tinuum must escape from galaxies at very high-z 
in order to explain the reionization of the universe 
which apparently was complete by z ~ 6 (Fan et 
al. 2006). 

Escaping Lyman continuum photons would not 
contribute to photoionization within galaxies and 
would thus lower the luminosities and equivalent 
widths of Lya and other Hydrogen recombination 
lines. However, Lyman continuum leakage would 
not necessarily lower the Lya escape fractions or 
the Lya/Ha ratios. On the contrary, the existence 
of Hi clear paths through the ISM would also pro- 



galaxies and their result suggest that in starbursts, 
i.e. galaxies with high Ha surface brightness, some 
of the Lyman continuum radiation might be leak- 
ing out, and there is tentative direct detection of 
Lyman continuum emission in Haroll (see Sect. 
4.5). At high rcdshift, determining direct escape 
fractions really means pushing the observational 
capabilities, and implied escape fractions are in 
general low (Siana et al. 20071, although several 



galaxies have been found with significant Lyman 
continuum leakage (Shapley et al. 2006). 



6. Perspectives 

Not only is Lya currently the most favorable 
window on the galaxy population at high-z but 
it is expected to remain so. Balmer lines are 
naturally much easier to interpret than Lya but 
their longer rest-wavelength makes them less suit- 
able for high-z. The James Webb Space Tele- 
scope promises huge improvements, particularly 
regarding Ha imaging and spectroscopy at high- 
z and near- and mid-IR continuum observations 
from an environment free from atmospheric back- 
ground emission. However, at z > 6, Ha will be 
redshifted out of the bandpasses of NIRCAM and 
NIRSPEC, leaving only MIRI with limited field- 
of-view and multiplexing capabilities for the de- 
tection of Ha. Thus especially if Lya is, as pre- 
dicted, visible from beyond the epoch of reioniza- 
tion (e.g. Cen & Haiman 2000) then Lya is likely 
to remain one of the primary probes of the high 
redshift universe for the coming decades. 

Some rather interesting possible trends have 
been identified and to improve on the study pre- 
sented here a number of alternatives can be en- 
visaged. Comparison of the co-evolution of H^ya 
and Wn a yields results that are difficult to attach 
meaning to in a small, hand-picked sample. What 
is required to rectify this is an imaging survey of 
nearby star-forming galaxies in Lya, selected in 
a clean homogeneous manner. Some of the tar- 
gets presented here are of sufficient star-formation 
rate and FUV luminosity that they may be di- 
rectly comparable to some of the objects being 
studied by high-z Lya observations and such a 
low-z sample would prove extremely valuable for 
co-interpretation of the survey data. Indeed, sim- 
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ulations for such a study have already been per- 
formed (Paper i) and have demonstrated that re- 
sults of similar reliability to those presented here 
may be obtained with a slightly reduced number 
of filters. 

The use of HST for this project was initially mo- 
tivated by the need for a FUV-capable space-based 
imager to observe Lya at z ~ and not primar- 
ily driven by the high spatial resolution. Given 
that we see variations in Lya between emission 
and absorption on scales down to the resolution 
limit of the telescope, the need for the high re- 
solving power has become essential. H/3 obtained 
with HST would permit a direct examination of 
extinction in the gas phase, an analysis of the role 
of dust on a pixel-by-pixel basis, a better evalua- 
tion of the ionizing photon production, and would 
also enable us to lock the nebular E(B — V) in the 
continuum fitting software. An alternative that we 
are currently exploring is the use of ground-based 
adaptive optics NIR imaging in the Br7 line using 
the NAOS-CONICA instrument on ESO's VLT. 

7. Conclusions 

We have produced continuum subtracted and 
calibrated Lya and Ha images for a sample of 
six local star-forming galaxies using the Advanced 
Camera for Surveys onboard Hubble Space Tele- 
scope. Morphological and photometric analyses 
have been performed on the Lya, Ha and stellar 
continuum maps. In summary: 

• A variety of morphologies is found, ranging 
from general absorption to emission. Ab- 
sorption and emission in individual targets 
are found to change on scales down to the 
resolution limit of the telescope. Lya fea- 
tures in emission and absorption show only 
minor correlation with the stellar morphol- 
ogy 

• Typically Lya emission is seen on large 
scales surrounding the central star forming 
regions in the form of low surface bright- 
ness halos. Moreover, Lya is found to be 
more spatially extended than Ha and ap- 
pears to emerge almost preferentially at low 
FUV surface brightnesses. At fainter isopho- 
tal levels, Lya is typically seen in emission, 
with fluxes and equivalent widths greater 



than those that would be predicted from 
recombination theory. 

• Observed Lya escape fractions range from 
below (absorber) to about 14%. Normal 
corrections for internal dust extinction are 
insufficient to explain these values. 

• Although any trends found from such a small 
hand picked sample must be viewed with 
caution, the found anti-correlation between 
the equivalent widths of Lya and Ha is quite 
suggestive. This could be explained if Lya 
escape is regulated by feedback, i.e. the de- 
velopment of wind blown bubbles and ISM 
flows that allow the photons to avoid the res- 
onant trapping by static Hi. 

• All the above statements point directly to 
one conclusion: the phenomenon of resonant 
scattering is observed to be extremely im- 
portant in the regulation of Lya emission, 
and radiative transfer of Lya photons in 
star-forming galaxies. 

• A few of our galaxies have far ultraviolet lu- 
minosities comparable to some of the objects 
being studied through their Lya emission at 
the highest redshifts. Given the low Lya es- 
cape fractions and complex resonant scat- 
tering phenomenon under observation here, 
care should be executed in interpretation of 
high-z survey data. 

• Our sample is fully calibrated and of physi- 
cal resolution three orders of magnitude bet- 
ter than the observed candidates at high-z. 
It may therefore be valuable in interpreta- 
tion of high-z datasets and of use to compare 
with the output of radiative transfer codes. 
The images are being released to the com- 
munity. 

As a service to the community we release cal- 
ibrated Lya, Ha, and their respective continuum 
images to the public domain, for the purpose of 
modeling high-z galaxies and Lya radiative trans- 
fer. The images and documentation are available 
at Centre de Donnees astronomiques de Stras- 
bourg (CDS) and also at the following website: 
http: //ttt . astro . su. se/projects/Lyman-alpha/ 
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Other data products, such as maps of age and red- 
dening and Lyman continuum production will be 
available upon request. 

GO and MH acknowledge the support from the 
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Research Council (Vetensapsradet) . We thank D. 
Schaerer for useful discussions. 
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